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Abstract 

In a helical flow there is a subrange of the inertial range in which there is a cascade 
of both energy and helicity. In this range the scaling exponents associated with 
the cascade of helicity can be defined. These scaling exponents are calculated from 
a simulation of the GOY shell model. The scaling exponents for even moments 
are associated with the scaling of the symmetric part of the probability density 
functions while the odd moments are associated with the anti-symmetric part of 
the probability density functions. 

Key words: anomalous scaling, helicity, cascade, shell models 
PACS: 47.27.Gs, 47.27.Jv 



1 Introduction 



In a helical flow both energy and helicity are inviscid invariants which are 
cascaded from the integral scale to the dissipation scale [1] . If these scales for 
the helicity are separated there will be an inertial range in which an equivalent 
of the four-fifth law for helicity transfer holds. This is a scaling relation for 
a third order structure function with a different tensorial structure from the 
structure function associated with the flux of energy. For helicity flux this is, 
(5v||(Z) • [v_|_(r) x v±(r + I)]) = (2/15)<5/ 2 , where 5 is the mean dissipation of 
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helicity. This relation is called the 'two-fifteenth law' due to the numerical 
prefactor [2,3]. The inertial ranges for helicity cascade and for energy cascade 
are different because the dissipation of helicity scales as D H {k) ~ kD E (k), thus 
the helicity will be dissipated within the inertial range for energy cascade. 
From balancing the helicity flux and the helicity dissipation a Kolmogorov 
scale £ = K H l for helicity dissipation can be defined [4], 

e ~ {^/fy\ (i) 



where v is the kinematic viscosity, e is the mean energy dissipation per unit 
mass and S is the mean helicity dissipation per unit mass. This scale is larger 
than the usual Kolmogorov scale rj = K^ 1 ~ (^ 3 /^) 
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The physical picture for fully developed helical turbulence is shown schemat- 
ically in figure 1. The mean dissipations 5 and e are solely determined by 
the forcing in the integral scale. There will then be an inertial range with 
coexisting cascades of energy and helicity with third order structure functions 
determined by the four-fifth - and the two-fifteenth laws. This is followed by 
an inertial range between Kh and Ke corresponding to non-helical turbu- 
lence, where the dissipation of positive and negative helicity vortices balance 
and the two-fifteenth law is not applicable. 
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Fig. 1. The inertial range for helicity cascade is smaller than the inertial range for 
energy cascade. In the range Kh < k < Ke the dissipation of positive and negative 
helicity balance. 
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2 The anomalous scaling exponents 

There is now experimental evidence that the K41 scaling relations are not 
exact. There are corrections for moments different from 3, expressed through 
anomalous scaling exponents, (5v(l) p ) ~ where 7^ p/3. Understand- 
ing and quantitatively determining the anomalous scaling exponents is one of 
the most intriguing and unsolved problems in turbulence. The intermittency 
corrections to the K41 scaling could depend on the transfer of helicity, maybe 
similar to the way the different sectors in anisotropic turbulence might give 
rise to sub-leading corrections of scaling exponents [5]. Furthermore, the he- 
licity cascade itself leads to a set of anomalous scaling exponents related to 
moments of the third order correlator of the two-fifteenth law. There is at 
present no experimental measurements from helical turbulence of the scaling 
exponents associated with the two-fifteenth law. 

It was recently shown numerically by Biferale et al. [6] that in the case of a 
shell model the anomalous scaling exponents for the helicity transfer has a 
strong difference between odd and even powers such that the scaling exponent 
( H (p) is not a convex function. 

Biferale et al. used a shell model consisting of two coupled GOY shell models. 
We will show here that the results obtained holds for the standard GOY shell 
model as well. Shell models are toy-models of turbulence which by construc- 
tion have second order inviscid invariants similar to energy and helicity in 3D 
turbulence. Shell models can be investigated numerically for high Reynolds 
numbers, in contrast to the Navier-Stokes equation, so that high order statis- 
tics and anomalous scaling exponents are easily accessible. Shell models lack 
any spatial structures so we stress that only certain aspects of the turbulent 
cascades have meaningful analogies in the shell models. This should especially 
be kept in mind when studying helicity which is intimately linked to spatial 
structures, and the dissipation of helicity to reconnection of vortex tubes [7]. 
The following thus only concerns the spectral aspects of the helicity and energy 
cascades. 

The GOY model [8-10] is the most well studied shell model. It is defined from 
the governing equation, 

u n = ik n (u n+2 u n+1 - ^M n+ iM n _i + t —^-u n _ x u n -'i)* - vk 2 n u n + /„ (2) 

with n = l,...,N where the the complex shell velocities. The wave 

numbers are defined as k n = A™, where A is the shell spacing. The second and 
third terms are dissipation and forcing. The model has two inviscid invariants, 
energy, E = Y m E n = T in Kl 2 > and 'helicity', H — YJ n H n — £ n (e- iy n \u n \ 2 . 
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The model has two free parameters, A and e. The 'helicity' only has the correct 
dimension of helicity if |e — l|~ n = k n =>- 1/(1 — e) = A. In this work we use 
the standard parameters (e, A) = (1/2, 2) for the GOY model. 

A natural way to define the structure functions of moment p is through the 
transfer rates of the inviscid invariants, 



S*(K) = {^ E n ) p,3 )K p/3 - K cE(p) (3) 
s; {K) = ((^) p/3 )K 2p/3 - K cH(p) (4) 

The energy flux is defined in the usual way as ITf = d/dt\ n ,i(J2m=i E m ) where 
djdt\ n x. is the time rate of change due to the non-linear term in (2). The 
helicity flux IT^ is defined similarly. By a simple algebra we have the following 
expression for the fluxes, 



(ITf ) = (1 - e)A n + A n+1 = S (5) 
(n^) = (-l) n k n (A n+1 - A n ) = S (6) 

where A n = k n _ilm{u n -iu n u n+ i) , e and 8 are the mean dissipations of energy 
and helicity respectively. The first equalities hold without averaging as well. 
These equations are the shell model equivalents of the four-fifth - and the 
two-fifteenth law. 

In the definition (3), (4) of the structure functions there is a slight ambiguity 
in the definition of x p ^ 3 for negative x and p not a multiplum of 3. The complex 
roots for (-1) 1/3 are (-1, 1/2 ± iy/3/2) and for (-1) 2 / 3 they are (1, -1/2 ± 
iy/3/2). The common way of circumventing the ambiguity is by defining x p * = 
sgn(x) |x| P//3 , which neglects the imaginary rootsQ. With this definition we 
have, 



S p {k n ) = J ty n (x) + M-x)V /3 dx = J *P+{x)x p / 3 dx (p even) (7) 



oo oo 

S P (k n )= l[^ n {x) -ip n (-x)]x p/3 dx = [ ip-(x)x p/3 dx (p odd) (8) 



where ip n {x) is the probability density function (pdf) for n^. ipni x ) ^ s (twice) 
the symmetric part of the pdf and ip~ (x) is (twice) the anti-symmetric part. 
Note that ipn( x ) * s itself a pdf while ipn( x ) * s n °t- i>n ( x ) exce Pt f° r a 
normalization, a pdf only if ip n {x) > ip n (—x) for all positive x. 



2 This can not always be done. Had we defined the structure functions from some, 
say, sixth order correlator, we are in trouble since z = (— l) 1 ' 6 has no real roots. 
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Fig. 2. The anomalous scaling exponents, C E (p) (dashed curve) and C, H {p)- The two 
full curves for even and odd moments for the helicity are calculated according to (7) 
and (8) respectively. The error bars on the ( E (p) curve are within the triangles, so 
we see a deviation between C, E {2p) and (, H (2p). The error bars for the odd moments 
of C H {p) are large because they are determined by the scaling of the anti-symmetric 
part of the probability density functions. The dashed-dotted line is the K41 scaling. 

The scaling exponents are determined from the scaling of the pdf's through, 

oo oo 

(x)dx, (9) 

— oo — oo 



so the scaling exponents for p even is related to the scaling of ip+ while for 
p odd they are related to the scaling of We have performed a simulation 
of the standard GOY model with (e, A, v, N) = (1/2, 2, 1CT 9 , 26) and a forcing 
of the form f n = 0.015 2 ,n/w2, corresponding to a constant energy input. The 
simulation was about 5000 large eddy turnover times. 

Figure 2 shows the anomalous scaling exponents, ( E {p),( H {p) f° r the energy 
and the helicity calculated according to (7) and (8). Using (7) the scaling ex- 
ponent ( H (p) can be denned for any real positive p, which from the Holder in- 
equality is a convex curve. Similarly using (8) assuming <f>~(x) to be a positive 
function we can define a continuous curve ( H (p) which is also from the Holder 
inequality a convex curve. The scaling exponent ( H (p) defined for integer p 
jumps between the two curves shown in figure 2. The scaling exponents differs 
from the ones found by Biferale et al. for the two-component GOY model. We 
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find that ( H (2p) is slightly larger than ( E (2p). The scaling regime in which 
( H (p) is calculated is < k < while for the energy K-[ < k < Ke- 
The negative part of the probability density is negligible in the case of energy 
transfer, ip n (x) ~ tp^ (x) ^ni x ) f° r ^ > 0, but for helicity transfer the neg- 
ative tail is big which gives the strong even-odd oscillations between the two 
curves. Note that C^(3) = 1 and ( H (3) = 2 are just the four-fifth - and the 
two-fifteenth law. 
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Fig. 3. The left hand panels show probability density functions ip(x) of the helic- 
ity flux for shells 3 and 6. The right hand panels show the probability distribution 
functions for x > 0, the dashed curves are 1 — x) which for a symmet- 

ric distribution is identical to ^(x). The scaling exponents for the even moments 
depends on the scaling with n of the symmetric part of the pdf which corresponds 
to the mean curve between the full and dashed curves. The scaling exponents for 
the odd moments depends on the scaling of the anti-symmetric part of the pdf 
corresponding to the gap betwwen the full and dashed curves. 

Figure 3 shows the probability distribution function (PDF) for helicity flux, 
defined by fy(x) = Jf ^ ip(y)dy for shell numbers n = 3 and n = 6 both in the 
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inertial range for helicity flux. The negative tail is plotted as 1 — x) which 
for a symmetric pdf gives two overlapping curves. We can similarly define the 
PDF *f?n(x) = Jq ip±(y)dy. A simple algebra gives = * n (x) + (1 - 

*„(-z))-l and *"(a;) = * n (z) - (l-* n (-x)) + (l -2*(0)). So we see that 
the scaling of ^^(x) is related to the scaling of the mean of the two curves in 
the right panels in figure 3, while the scaling of ^(x) is related to the gap 
between the two curves. 



3 Conclusion 



Coexisting cascades of energy and helicity are possible in the GOY shell model. 
The scaling of the odd order moments of the helicity transfer depends on the 
scaling of the anti-symmetric part of the probability density function for the 
helicity flux. This defines a convex anomalous scaling curve through the point 
C (3) = 2 which is the two-fifteenth law. The even order moments of the 
helicity flux has anomalous scaling exponents close to the ones found for the 
energy flux. In the simulation a scale break at Kh is not observed. This implies 
that the anomalous scaling exponents for the energy flux are not influenced 
by the cascade of helicity. 
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